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Chemical functionalization of graphene holds promise for various applications ranging from nano- 
electronics to catalysis, drug delivery, and nanoassembly. In many applications it is important to be 
able to transport adsorbates on graphene in real time. We propose to use electromigration to drive 
the adsorbate transport across the graphene sheet. To assess the efficiency of electromigration, we 
develop a tight-binding model of electromigration of an adsorbate on graphene and obtain simple 
analytical expressions for different contributions to the electromigration force. Using experimen- 
tally accessible parameters of realistic graphene-based devices as well as electronic structure theory 
calculations to parametrize the developed model, we argue that electromigration on graphene can 
be efficient. As an example, we show that the drift velocity of atomic oxygen covalently bound to 
graphene can reach ~1 cm/s. 



Many unique properties of graphene - a monoatomic 
crystalline sheet of carbon - stem from the fact that it 
"lacks" volume and is, therefore, a truly 2d "all-surface" 
material pQ. For instance, the surface functionalization 
of graphene (e.g., graphene oxide) provides an oppor- 
tunity to alter electronic properties of the entire mate- 
rial, which holds promise in nanoelectronics [HE], non- 
volatile memory [4], graphene-based nanoassemblies for 
catalysis, photovoltaics and fuel cells applications [5j|6]. 
In majority of these applications, the performance of a 
graphene-based device can be significantly improved pro- 
vided there is a way to tune the surface functionalization 
in real time, i.e., during device operation. The related 
problem is to control and direct the transport of adsorbed 
atoms/molecules for nanoassembly [7 and drug delivery 
applications [8]. 

Electromigration is the drift of material on the surface 
(or in the bulk) of a current-carrying conductor [9 . We 
propose to exploit electromigration as an efficient and 
easily controllable method to drive the directed trans- 
port of adsorbates on graphene. Very recently, and for 
the first time, the efficient electromigration of metallic 
clusters/atoms on graphene has been demonstrated ex- 
perimentally [10]. In the current work, however, we fo- 
cus on a different class of adsorbates (perhaps more rele- 
vant from the perspective of chemical functionalization) 
- atoms or molecules covalently bound to graphene. At 
first glance, the covalent binding may seem too strong to 
allow for efficient electromigration. Nevertheless, recent 
theoretical studies suggest the possibility of fast diffu- 
sion. Specifically, the activation energy for the atomic 
oxygen diffusion has been shown to drop from ~0.7 eV 
for neutral graphene to ~0.15 eV for n-doped graphene, 
resulting in a diffusion coefficient as high as ~ 10 -6 cm 2 /s 

HB. 

In this Letter, we investigate electromigration of an ad- 
sorbate covalently bound to graphene, focusing on atomic 
oxygen (O) and the amino group (NH), as an exam- 
ple. In the lowest energy configuration, a single oxygen 




FIG. 1. Electromigration on graphene. (a) Schematic de- 
piction of the electron electromigration force on an adsorbate 
(oxygen atom) bound to graphene surface: direct force Fd 
and the force due to an electron scattering event, £F k / )k , 
contributing to electron wind force F^. The adsorbate is 
negatively charged and both Fd and F^ are directed oppo- 
site to the current. The current flows downward, as shown. 
(b),(c),(d) Adsorbate hopping: the adsorbate remains bound 
while hopping (e.g., from position 1-2 to position 2-3). The 
arrow indicates the net drift direction. The highlighted bonds 
in inset (d) are excluded from the tight-binding description 
due to the presence of a covalently bound adsorbate (see text) . 



atom covalently binds to two adjacent carbon atoms of 
the graphene's honeycomb lattice, thus forming an epoxy 
bond, Figs. [[{a^d). 

A single nitrogen atom can bind similarly, providing its 
third covalent bond for chemical functionalization, e.g., 
in drug delivery applications. In what follows, we develop 
a simple tight-binding model and obtain an analytical re- 
sult for the drift velocity of an adsorbate as a function 
of its charge, the electric current in graphene, as well 
as the backgate doping level, and temperature. Specifi- 
cally, we find that the migration (drift) velocity reaches 
up to ~ 0.6 — 4 cm/s at electrical current densities of 
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FIG. 2. Diagrams representing various contributions to the 
electromigration force (see text for details). Zigzag lines de- 
note the external dc electric field, crosses stand for the ad- 
sorbate scattering vertexes. Dashed lines are the adsorbate- 
substrate interaction. The grey-filled blob denotes the elec- 
tron density in graphene. Solid lines stand for the equilib- 
rium propagator of electrons in graphene. Thick solid lines 
denote fully dressed nonequilibrium electron propagators (in 
the presence of the electric field). Double dashed lines repre- 
sent the screened adsorbate-substrate interaction. 

~ 1 A/mm and temperatures of 300 — 500 K. The doping 
level and temperature dependence of the drift velocity are 
then used to formulate a robust adsorbate manipulation 
technique based on local heating of a graphene sample. 
This technique can become suitable for patterning and 
other applications involving direct access to the surface 
of graphene. 

A particle (or defect) in contact with a conductor ex- 
periences an electromigration force which can be written 

as [9] 

F = F d + F w =eZE-<V R £7), (1) 

where e, Z and E are the absolute value of the electron 
charge, the charge of the particle in atomic units, and the 
vector of the external dc electric field (EF), respectively. 
The first term, the direct force F^, originates from the 
direct interaction of the adsorbate's charge with the EF, 
as shown diagrammatically in Fig. [2|a). 

The second term, F w = -(VrU) [Fig. |^b)], orig- 
inates from the scattering of electrons in graphene by 
the adsorbate-graphene interaction potential, U (R is the 
position of the adsorbate). This contribution is often re- 
ferred to as the electron wind force. Assuming the adsor- 
bate structureless, i.e, treating U as just a single-particle 
scattering potential, this contribution can be rewritten 
as 

F w = «Fk',k =\Y. fi ( k ' - k )£Wkk', (2) 

k,k' k,k' 

where pkk' = (c^/Ck) is the single-particle density ma- 
trix and c£ (£k) creates (annihilates) a Bloch wave with 
quasimomentum k. Equation ^ has a very appealing 
microscopic interpretation: it describes electron scatter- 
ing from k-state to k'-state within graphene, accompa- 
nied by fi(k — k') momentum transfer to the adsorbate, 
with amplitude C/k'kPkk', see Fig.[TJa). The pitfall is that 
Eq. [2] is exact only for free electrons, and it can be non- 
trivial to justify that the momentum transfer is equal to 



the change of quasimomentum in the presence of band- 
structure effects [9]. In what follows, we will show that 
the dominant contribution to U is sufficiently smooth (as 
compared to the size of the graphene unit cell), so that 
Eq. ([2} holds. 

To obtain a closed expression for F w , the density ma- 
trix pkk' in Eq. Q is expanded up to the leading order 
in U. The validity of such perturbative expansion will be 
justified later when we address the graphene-specific form 
of U. The first-order contribution in U to F w vanishes 
exactly. The simplest nonvanishing contribution (second 
order in U, first order in E) describes the interaction of 
an adsorbate with the current-carrying charge density of 
graphene, Fig. [2|c). Important processes not described 
by this diagram include (i) Coulomb-induced screening 
of U and (ii) scattering due to impurities and phonons. 
Most diagrams containing such processes, as well as the 
effect of the EF, can be "lumped together" by dress- 
ing the adsorbate-graphene interaction using the random 
phase approximation and replacing bare electron propa- 
gators with fully dressed nonequilibrium (i.e., current- 
carrying) Green's functions. As a result, the electron 
wind contribution to the electromigration force, depicted 
diagrammatically in Fig. |2jd), can be expressed as 

F w = -^h(k!-k)ST^ (3) 

k,k' 

<5IV k = F |C/ k , k | 2 [l - / i (k')]/ J (k)5( Wfe , fc ), (4) 

where L^'k is the screened adsorbate-graphene poten- 
tial. The steady-state distribution function in graphene 
is given by 

f j ( k ) = f (k)- 4 ^^5(k-k F ), (5) 

in the linear response approximation [12] . Here j is the 
current density within graphene and fo(k) is the Fermi- 
Dirac distribution; the Fermi momentum and the velocity 
of Dirac electrons in graphene are denoted by kp and 
vf, respectively. Equations Q and Q agree with the 
general (second-order) result for the wind force obtained 
earlier in Ref. [9 j. To complete the derivations we need 
to obtain /7kk', which is specific to the case of adsorbates 
on graphene. 

The interaction potential between an adsorbate and 
graphene has two distinct contributions: U = U c + U dei . 
The first one, U c , describes the scattering of electrons in 
graphene by the Coulomb potential of the charged adsor- 
bate. The second one, £7 def , describes the scattering due 
to the lattice defect caused by the covalent bonding of 
an adsorbate to graphene. As will become clear shortly, 
the Coulomb contribution is dominant and, therefore, we 
discuss it first. 
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In the backgated graphene with experimentally acces- 
sible electron densities, the screened Coulomb poten- 
tial varies slowly over graphene unit cell, and it suf- 
fices to consider scattering only within a single Dirac 
cone. Within the lowest energy tight-binding descrip- 
tion and the random phase approximation for screening, 
the Coulomb contribution can be cast in the form [13] [T4] 

2ire 2 Z/k (1 + 



^k'k 



k 



-|k'-k|/i 



<c TF 



(6) 



where krF = ^e 2 k F / '(hkv F ) ~ 9k F /h is the Thomas- 
Fermi momentum. For graphene laid on top of a half- 
space dielectric substrate with the dielectric constant ft, 
the effective constant is given by R = (ft + l)/2 [T5] . 
In what follows, a Si02 substrate with ft = 2.5 will be 
assumed. The angle ^kk' stands for the angle between 
the vectors k and k', both measured with respect to the 
same Dirac point. The distance between the adsorbate 
and graphene, h, is comparable to the graphene lattice 
constant, a, and, hence, |k' — k.\h< 2k F h <C 1 is negligi- 
ble. 

The effective size of the potential U c is given by 
fc^ F , so that at \i = 0.2 eV (easily reached in back- 
gated graphene) one has k F a ~ 0.08, and, therefore, 
krF(i ~ 0.3. Therefore, the size of the potential is signif- 
icantly larger than the unit cell, justifying both Eq. (|3| 
and the neglect of scattering between the Dirac cones. 

Using the inequality krF ^ fcp> the screened Coulomb 
potential of the adsorbate simplifies to /7^ k « ne 2 Z(l + 
e 26>kk/ )/ft/cTF- Substituting this expression, along with 
Eq. ([5]), into Eq. Q and retaining only terms linear in 
j, we finally obtain for the Coulomb contribution to the 
electron wind force [16] 



Hj ( 7rZ 



(7) 



Here, the spin and valley degeneracies are lumped into 
the current j. 

To prove that the Coulomb contribution to the elec- 
tron wind force is dominant, we estimate the amplitudes 
of Coulomb and lattice defect contributions. The co- 
valent binding of an adsorbate to graphene alters the 
hybridization of the involved carbon atoms from sp 2 to 
sp 3 , see Figs.[TJb)-[TJd). This amounts to cutting out the 
corresponding p z carbon orbitals from the tight-binding 
description of graphene, creating a lattice point defect. 
The lack of a single p z orbital can be mimicked by a fic- 
titious impurity potential which cancels out the hopping 
integrals involving this orbital. For electrons near the 
Fermi circle such an impurity potential is essentially a 
delta function with a magnitude in momentum space of 
Z7 def ~ ta 2 ~ hvpCL, where t « 2.8 eV is the hopping en- 
ergy. Comparing U c and U def (screened or unscreened), 
we obtain 



U det /U c (k F ) 



hvj- 



2tt\Z 



k F CL : 



0.07 



Since k F a <C 1 (see above), the contribution of the single- 
bond defect potential is small compared to that of the 
Coulomb potential at not too small Z. The actual defect 
potential leads to exclusion of several bonds, e.g., five for 
oxygen in the equilibrium state; see Fig.[jjd). Further- 
more, the defect potential cannot be considered smooth 
on the scale of the unit cell, which might introduce signif- 
icant band-structure effects [see the discussion following 
Eq. Q]. Nevertheless, our analysis (not provided) using 
the accurate tight-binding model showed that these two 
complications do not change the qualitative conclusion 
that we have drawn: the U def contribution to the elec- 
tron wind force is small compared to the Coulomb one 
and will be omitted henceforth. 

Finally we note that the second-order representation of 
the wind force, Eq. Q and Fig.[2|d), is sufficient for ad- 
sorbates with not very large Z. As will be seen shortly, 
\Z\ <0.4 for both oxygen and nitrogen, and, therefore, 
the screened Coulomb potential is sufficiently small rela- 
tive to the Fermi energy U c (k F )k 2 F /v F k F « \Z\/2 < 0.2. 
At larger Z, multiple events of adsorbate-induced scat- 
tering of in-graphene electrons can modify the results ob- 
tained here. 

The total driving force of the electromigration for a 
charged adsorbate on graphene is given by 



Hj ( 7rZ 



eZ*j/a, (9) 



k F a. (8) 



where the effective charge of the adsorbate is introduced 

as Z* = Z — and <7n = 1% is the minimal conduc- 

1 — 1 
tivity of graphene [14] . Equation rt9l) - the main result of 

this work - expresses the total driving force of the electro- 
migration via experimentally accessible parameters such 
as conductivity of graphene and the current density. The 
charge of the adsorbate, however, is not directly accessi- 
ble, so we estimated it by performing electronic structure 
theory calculations using Gaussian 09 quantum chem- 
istry package [17 . These calculations were performed 
using density functional theory with the PW91 functional 
[T8] and the 6-31 1G** basis set for a number of graphene 
flakes of increasing size (up to C62) to guarantee the con- 
vergence with respect to boundary effects. A single ad- 
sorbate was put in the center of a flake. The charges of 
the adsorbates were found in the range Zo G (—0.2,-0.4) 
and Znh £ (—0.1,-0.3), depending on the specific po- 
sition of the adsorbate along the hopping trajectory [see 
Fig.[TJb)-[TJd)], and the graphene doping level. The ratio 
of the electron wind and direct forces (or, equivalent ly, 
Z*/Z — l) is shown in Fig.|3jb) (as a function of graphene 
conductivity [19J) for the range of Z values. As is seen, 
the electron wind contribution dominates the electromi- 
gration force except for smallest Z or low conductivity,. 

To calculate the drift velocity of an adsorbate on 
graphene due to electromigration, we use the Einstein- 
Smoluchowski relation between the diffusion coefficient, 
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FIG. 3. Drift of an adsorbate along the graphene sheet, 
(a) Drift velocity, v, for Z=-0.3 j = 1 A/mm, a — 1 mO" 1 
and different temperatures. Dashed vertical lines mark acti- 
vation energies of 0.15 eV and 0.7 eV, which correspond to 
oxygen diffusion on n-doped and charge- neutral graphene, re- 
spectively [11] . (b) Relative contribution of the electron wind 
force to the effective charge of the adsorbate. (c) Robust con- 
trol of activation and drift path with a local heater, e.g. laser 
beam. 

and the drift velocity, v, i.e., 

v = FD/k B T, D = ^° e -JW*BT (1Q) 

where d = 1.23 A is the hopping distance. The attempt 
frequency for oxygen was found to be v$ = 26 THz in 
Ref. [11 . We assume the same attempt frequency for 
the NH group, having in mind the approximate char- 
acter of the calculations and the fact that the diffu- 
sion coefficient is not overly sensitive to variations of 
this parameter (compared to, e.g., the activation en- 
ergy or temperature). The activation energy for oxy- 
gen diffusion, E ai as obtained in Ref. [11], is ~0.7 eV 
for charge- neutral graphene, and a much lower value of 
~0.15 eV for n-doped graphene with the charge density 
of —7.6 x 10 13 cm -2 . Our electronic structure theory 
calculations qualitatively confirm the strong sensitivity 
of the activation energy for oxygen diffusion. Further- 
more, we found a similar dependence of the activation 
energy on doping level for NH adsorbate, albeit with 
somewhat higher activation energies (by 0.2 — 0.3 eV). 
Figure [3^ a) shows the dependence of the drift velocity on 
E a for Z=-0.3 at j=lA/mm and a=10" 3 Or 1 . Specifi- 
cally, at these parameters and £? a =0.15 eV the drift ve- 
locity is ~6 mm/s at room temperature (T=300 K) and 
reaches up to 4 cm/s at T=500 K. 

The increased drift velocity of the electromigration at 
higher temperatures can be used to perform adsorbate 
manipulations and patterning via a guided motion of ad- 
sorbates along the surface of graphene. This is most 
easily achieved by heating graphene locally with a fo- 
cused laser beam [Fig. [3^c)], or by a heated AFM tip 
[20] . This local heating will enhance u, while keeping the 



drift velocity low outside the heating spot by adjusting 
the doping level to lower values. As a result, it should 
be possible to move adsorbates along the desired path 
and assemble them into desired patterns by tracing the 
motion with the "local heater" and choosing the appro- 
priate current directions via a set of source-drain con- 
tacts in perpendicular directions (typical, e.g., for Hall 
conductance measurements), see Fig^c). The spatial 
resolution of patterning is thermodynamically limited by 
the balance between diffusion and drift, and is given by 
l p = kT/F [21]. For the current density and conductiv- 
ity used in Fig ^a) we obtain l p = 0.09T/Z* nm K _1 . 
For adsorbed oxygen atoms at T = 300 K, the result in 
l p ~ 50 nm. 

An experimental verification of an efficient electromi- 
gration of adsorbates on graphene can either be done 
using AFM/STM techniques to directly monitor diffu- 
sion/drift of adsorbates, or by optical means, e.g., adopt- 
ing a nitrogen-based adsorbate with a fluorescent func- 
tional group and monitoring the fluorescence of such ad- 
sorbates with temporal and spatial resolution. Mapping 
a trajectory of adsorbates on carbon materials with finite 
band gap (e.g., carbon nanotubes) by photoluminescence 
quenching is an alternative strategy [22] . 
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